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ABSTRACT 

The outer shells of young supernova remnants (SNRs) are the most plausible 
acceleration sites of high-energy electrons with the diffusive shock acceleration 
(DSA) mechanism. We studied spatial and spectral properties close to the shock 
fronts in four historical SNRs (Cas A, Kepler's remnant, Tycho's remnant, and 
RCW 86) with excellent spatial resolution of Chandra. In all of the SNRs, hard 
X-ray emissions were found on the rims of the SNRs, which concentrate in very 
narrow regions (so-called "filaments"); apparent scale widths on the upstream 
side are below or in the order of the point spread function of Chandra, while 
0.5-40 arcsec (0.01-0.4 pc) on the downstream side with most reliable distances. 
The spectra of these filaments can be fitted with both thermal and nonthermal 
(power-law and SRCUT) models. The former requires unrealistic high temperature 
(>2 keV) and low abundances (<1 solar) for emission from young SNRs and may 
be thus unlikely. The latter reproduces the spectra with best-fit photon indices 
of 2.1-3.8, or roll-off frequencies of (0.1-28) x 10^^ Hz, which reminds us of the 
synchrotron emission from electrons accelerated via DSA. We consider various 
physical parameters as functions of the SNR age, including the previous results 
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on SN 1006 (Bamba et al. 2003b); the filament width on the downstream side 
increases with the SNR age, and the spectrum becomes softer keeping a non- 
thermal feature. It was also found that a function, that is the roll-off frequency 
divided by the square of the scale width on the downstream side, shows negative 
correlation with the age, which might provide us some information on the DSA 
theory. 

Subject headings: acceleration of particles — supernova remnants: individual 
(Cas A, Kepler's remnant, lycho's remnant, SN 1006, RCW 86) — X-rays: ISM 

1. Introduction 

Ever since the discovery of cosmic rays (Hess 1912), the origin and the acceleration 
mechanism up to more than ~TeV have been long-standing problems. Koyama et al. (1995) 
discovered synchrotron X-rays from shock fronts of a supernova remnant (SNR), SN 1006, 
and suggested that shocks of SNRs are cosmic ray accelerators. Several SNRs have recently 
been categorized into synchrotron X-ray emitters (G347.3— 0.5: Koyama et al. (1997); Slane 
et al. (1999), RCW 86: Bamba, Tomida, & Koyama (2000); Borkowski et al. (2001b), and 
G266.6— 1.2: Slane et al. (2001)). Recently, a significant number of SNRs with synchrotron 
X-rays have been discovered by the ASCA Galactic plane survey (Bamba et al. 2001; Ueno 
et al. 2003; Bamba et al. 2003a; Yamaguchi et al. 2004). These discoveries provide good 
evidence for the cosmic ray acceleration at the shocked shell of SNRs. The most plausible 
process of the cosmic ray acceleration is the diffusive shock acceleration (DSA) (e.g. Bell 
1978; Blandford & Ostriker 1978; Drury 1983; Blandford & Eichler 1987; Jones & EUison 
1991; Malkov & Drury 2001), which can accelerate particles on the shock into a power- law 
distribution, similar to the observed spectrum of cosmic rays showering on the earth. 

Apart from the global success of DSA, there are still many remaining problems. We have 
not yet fully understood detailed but important information, such as the maximum energy 
of particles, the configuration of magnetic fields, the injection efficiency from thermal plasma 
to accelerated particles, and the acceleration and de-acceleration history of particles around 
the shock fronts, and so on. This partly comes from the fact that there is no information 
about the spatial distribution of accelerated electrons of practical objects, which strongly 
reflects the above uncertainties. 

Bamba et al. (2003b) considered the spatial distribution of emission from accelerated 
particles as new information in order to understand the acceleration on the shock in SNRs; 
they observed the synchrotron X-ray emitting shell of SN 1006 with the excellent spatial 
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resolution of the Chandra X-ray observatory (see also Long et al. 2003) , and found that the 
emitting regions of synchrotron X-rays are incredibly thin ("filaments"), about 0.01-0.1 pc 
on the upstream side and 0.06-0.4 pc on the downstream side, at 2.18 kpc distance (Winkler 
et al. 2003). This result constrains the DSA theory and magnetic-field configurations; both 
an efficient acceleration scenario with strong magnetic field parallel to the shock normal and 
an inefficient one with a weak and perpendicular magnetic field are allowed (Yamazaki et al. 
2004b; Berezhko et al. 2003). 

In the present work, we searched for nonthermal filaments in other historical SNRs 
(Cas A, Kepler's remnant, Tycho's remnant, and RCW 86) and measured their scale widths 
and spectral parameters. It was also investigated whether or not the age of the SNR is 
correlated with the scale width and/or spectral parameters of filaments. If it is, such cor- 
relations might describe the time evolution of these parameters in a single SNR. Although 
the characteristics of the SNR should be considered, those correlations may also become 
clues to understand the time evolution of the maximum energy of electrons, the magnetic 
field configuration, the injection efficiency, and so on. This paper is organized as follows. 
We summarize the observation details of five SNRs with Chandra in §2. §3 briefly describes 
the results of the observations of the four SNRs. In §4 we discuss the origin of emission 
(§4.1), our analysis of individual SNRs (§4.2), the time evolution of the spatial and spectral 
parameters of filaments together with previous analysis of SN 1006, which have already been 
analyzed in Bamba et al. (2003b) (§4.3). We also briefly deal with the time evolution of 
energy densities of magnetic field and cosmic rays in this section. A summary of our results 
is given in §5. 



2 . O bser vat ions 

We used the Chandra archival data of the ACIS of five historical SNRs (Cas A, Kepler, 
Tycho, SN 1006, and RCW 86), as fisted in Table 1. Note that this data set includes almost 
historical shell-type SNRs. The satelfite and the instrument are described by Weisskopf et 
al. (2002) and Garmire et al. (2000), respectively. Data acquisition from the ACIS was made 

in the Timed- Exposure Faint mode. The data reductions and analysis were made using the 
Chandra Interactive Analysis of Observations (CIAO) software version 2.3. Using the Level 
2 processed events provided by the pipeline processing at the Chandra X-ray Center, we 
selected ASCA grades 0, 2, 3, 4, and 6, as the X-ray events. The effective exposure for each 
observation is also given in Table 1. 
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3. 



Individual Results 



3.1. Cas A (Gill. 7-2.1) 



Although the remnant's precise age is still uncertain, the most probable record of the 
explosion is SN 1680, made by Flamsteed (1725), which we adopt in this paper. The distance 
to the remnant was estimated to be 3.41°;^ kpc by Reed et al. (1995) using the measurement 
of optical proper motions. The average radio index of the whole remnant at 1 GHz is 0.78 
(Green 2004). 

Cas A emits hard X-rays up to greater than 10 keV, which were detected by RXTE (Allen 
et al. 1997), BeppoSAX (Favata et al. 1997), and OSSE onboard GRO (The et al. 1996). 
Vink & Laming (2003) has already found and analyzed one of the nonthermal filamentary 
structures with Chandra. There are many other filaments thinner than that which they 
analyzed. These filaments are considered to be profiles of a thin sheet-like structure, as 
suggested by Hester (1987); then, the scale width of these filaments is the upper-limit width 
of the real thickness in the radial direction of the sheets. Therefore, we consider other 
filaments in this paper, which are located much nearer to the aim-point (less than 25 arcsec) 
and are much thinner than those analyzed by Vink & Laming (2003). 

The upper- left panel of Figure 1 shows the south-eastern side of Cas A in the 5.0- 
10.0 keV band, in which there are a lot of filamentary structures. These filaments are more 
clumpy than those in SN 1006 (Bamba et al. 2003b). We selected two filaments as shown in 
the panel, which are straight enough to be made profiles and near the aim-point. 

Figures 2(a) and (b) show the profiles of two filaments. In order to estimate the scale 
width of these filaments, we used a simple, empirical fitting model for the apparent profiles, 
which is the same as used in the previous analysis in the SN 1006 case by Bamba et al. 



where A and xq arc the flux and position at the emission peak, respectively. The quantities 
Wu and Wd represent the e-folding widths on the upstream and the downstream sides of 
the apparent emission peak, respectively (hereinafter, "m" and "d" represent upstream and 
downstream sides, respectively). We treat the best- fit values smaller than 0.8 arcsec as the 
upper limit, because they highly suffer the infiuence of point spread function (PSF), which 
is about 0.5 arcsec within ~4 arcmin radius from the aim point. The fittings were made 
with profiles in the 5.0-10.0 keV, and were accepted statistically with the best-fit models 
and parameters shown in Figures 2(a) and (b), and Table 2. 



(2003b): 




on the upstream side 



on the downstream side. 



(1) 
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We then made the spectra of the filaments within the scale widths {xq — Wu < x < 
Xq + Wd) in Figures 2(a) and (b). The background spectra were made from the just off- 
filament downstream sides with the same width of the filaments. We compared the Si line 
photon counts with those in the continuum (5.0-10.0 keV) band in these regions, and found 
that photons in this band include thermal photons in a ratio of only less than 25%, under the 
assumption that all photons in the background regions are thermal, and both regions have a 
thermal component with the same physical parameters, except for their flux. Therefore, we 
treated all the 5.0-10.0 keV photons in our analysis of Cas A, as a part of the hard X-rays 
reported in previous observations (Vink & Laming 2003). 

The background-subtracted spectrum accumulated from the two filaments is shown in 
the upper-left panel of Figure 3, which is featureless with no line-like structure, and extends 
to the hard X-ray side, which was fitted with both an absorbed thin thermal plasma model in 
non-equilibrium ionization (NEI) calculated by Borkowski et al. (2001a) ("NEI" model), and 
an absorbed power-law model. The absorption column was subsequently calculated using 
the cross sections by Morrison & McCammon (1983) with the solar abundances (Anders & 
Grevesse 1989). Both models were well fitted statistically with best-fit curves and parame- 
ters, shown in Figure 3 and Table 3. We also applied the SRCUT model, which represents the 
synchrotron emission from electrons of power-law distribution with exponential roll-off in a 
homogeneous magnetic field (Reynolds & Keohane 1999). The spectral index at 1 GHz was 
fixed to be 0.78 according to a report by Green (2004). This model also well reproduces the 
spectra with best-fit values shown in Table 3. The spectra are much harder than those of 
the filament analyzed by Vink & Laming (2003). 



3.2. Kepler's Remnant (G4.5-|-6.8) 

Kepler's remnant is the second youngest SNR in our Galaxy recorded by human beings, 
which appeared in 1604 (Kepler 1606). DeLaney et al. (2002) observed this SNR with VLA, 
and found the spatial variation of spectral index from 0.85 to 0.6 with a mean value of 
0.7. The distance to the remnant is 4.8 ± 1.4 kpc using the H I observations with VLA 
(Reynoso & Goss 1999). Petre, Allen, & Hwang (1999) found hard X-rays (above 10 keV) 
from Kepler with RXTE, and XMM-Newton pointed out that at least the south-eastern edge 
emits nonthermal X-rays (Cassam-Chenai et al. 2003). 

The upper-right panel of Figure 1 shows a 4.0-10.0 keV band Chandra image of the 
eastern rim of Kepler. We can sec very sharp structures on the outer edge of the remnant, 
especially bright on the eastern side, which is similar to the filaments in SN 1006. In order 
to carry out spatial and spectral analysis of these structures, we selected two filaments, as 
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shown in the panel, which are <160 arcsec distant from the aim-point. Figures 2 (c) and (d) 
show the profiles of the filaments in the 4.0-10.0 keV band. They have sharp edges both on 
the upstream and downstream sides of the emission peak, and eq.(l) well represents their 
profiles with best-fit curves and parameters summarized in the panels and Table 2. 

We obtained the spectra of these filaments in the same way of the Cas A case. Back- 
ground regions were selected just from the downstream sides of the filaments, which are free 
from any other structures. The intensities of iron and silicon lines and continuum photons 
in the 4.0-10.0 keV band were compared, and it was found that more than 90% of the 
photons in the 4.0-10.0 keV band arc the hard X-ray component reported in previous obser- 
vations (Petre et al. 1999; Cassam-Chenai et al. 2003), with the same assumptions as in the 
Cas A case. Thus the derived scale widths represent those of the hard X-ray emission. The 
background-subtracted spectrum accumulated from these filaments is shown in the upper- 
right panel of Figure 3, which shows a hard and line-less feature. We fitted each spectrum 
with NEI and power-law models with an absorption. Both models are well fitted to the data 
with the best-fit models and parameters shown in the panel and Table 3. The SRCUT model 
was also applied with the spectral index at 1 GHz fixed to be 0.70, derived by DeLaney et 
al. (2002) as the average value on the eastern side. This model also well represents the data 
with the best-fit parameters listed in Table 3. 



3.3. Tycho's Remnant (G120.1-|-1.4) 

Tycho is a remnant of a supernova that exploded in 1572 (Tycho 1573). The distance to 
the remnant is estimated to be 1.5-3.1 kpc, from the proper motion of the optical filaments 
(Kamper & van den Bergh 1978) and the shock velocity estimated by Ghavamian et al. 
(2001). We adopt the most popular value, 2.3 kpc. Katz-Stone et al. (2000) investigated 
the variations of the radio spectral indices, and found that the emission in the outer rim 
shows a trend that brighter clumps have a flatter spectral index than the average index of 
s — 0.52 ± 0.02, maybe due to either SNR blast waves and ambient medium interactions or 
internal inhomogeneities of the magnetic field within the remnant. 

Hard X-ray emission has been unmistakably detected from Tycho, up to 30 keV, with 
HEAO 1 (Pravdo & Smith 1979), Ginga (Fink et al. 1994), and RXTE (Petre et al. 1999). 
Chandra and XMM-Newton also observed Tycho's remnant, and the good spatial resolutions 
of instruments onboard these satellites show the spatial variety of X-ray emission (Hwang et 
al. 2002; Decourchelle et al. 2001). Hwang et al. (2002) claimed that there is a thin structure 
encircling the SNR, the spectrum of which has lines with low equivalent widths. 
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In order to be clear as to whether the thin structure is a hard X-ray emitter or not, 
spatial and spectral analysis were conducted in the same way as in the cases for Cas A and 
Kepler. The lower-left panel of Figure 1 shows the north-western part of the remnant in 
the 2.0-10.0 keV band. The thin filament-like structure suggested by Hwang et al. (2002) 
can be seen. Figure 2 (e)-(i) shows profiles of the filaments selected in the panel, which 
are located within 120 arcsec from the aim-point. Very sharp rises and decays can be seen 
in all of the filaments. The fittings with eq.(l) were performed only around the outermost 
peak of hard X-ray emission, because thermal photons dominate in the downstream regions 
(the right side of each panel), and were accepted statistically with the best-fit models and 
parameters shown in these panels and Table 2. 

Their spectra were obtained in the same way as in the Cas A case. Background spectra 
were accumulated from just the downstream side of the filaments. The lower-left panel 
of Figure 3 shows the background-subtracted spectra, which were combined with those of 
filaments in each CCD chip. Since both are hard and have no line-like structure, we fitted 
them with NEI, power-law, and SRCUT models with an absorption. We fixed the spectral 
index at 1 GHz to be 0.52, following Katz-Stonc ct al. (2000). All three models were accepted 
statistically with similar probability, with the best-fit values listed in Table 3. We compared 
the Si and S lines and the 2.0-10.0 keV continuum band intensity in the source spectrum 
with those in background region, and it was found that only less than 6% of the thermal 
photons contribute in the 2.0-10.0 keV band. Therefore, it is hkely that the all photons in 
the 2.0-10.0 keV band can be treated as pure nonthermal. 



3.4. RCW 86 (G315.4-2.3) 

RCW 86 was identified as being the oldest historical SNR by Clark & Stephenson 
(1977). They suggested that a progenitor of the remnant was reported in AD 185 in Later 
Han dynasty records (Fan You 432), although Chin & Huang (1994) claimed some doubt 
about the record. The spectral index at 1 GHz and the distance was measured to be 0.6 by 
Caswell et al. (1975) and 2.8 kpc (Rosado et al. 1996), respectively. 

Strong hard X-rays have been detected from the remnant by Ginga (Kaastra et al. 
1992) and RXTE (Petre et al. 1999). Using ASCA data, Bamba et al. (2000) and Borkowski 
et al. (2001b) independently found that the hard X-rays show a nonthermal feature, and 
concentrate on the south-western (SW) shell of the SNR. The Chandra observation with 
excellent spatial resolution reveals that the regions emitting nonthermal X-rays are very 
clumpy and different from those emitting thermal X-rays (Rho et al. 2002). 
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The lower-right panel of Figure 1 represents a close-up view of the SW shell in the 
2.0-10.0 keV band. Clumpy structures can be seen, which have been already indicated as 
nonthermal filaments by Rho et al. (2002). We selected two filaments for our analysis, which 
are bright in the hard band image and are located <120 arcsec from the aim-point. Their 
profiles from the region indicated in the panel were made, as shown in Figures 2(j) and (k). 
They show a sharp rise and a rather slow decay in the upstream and the downstream sides of 
the emission peak, respectively. The fitting with eq.(l) was carried out in the 2.0-10.0 keV 
band, and was accepted statistically with best-fit models and the values shown in the panels 
and Table 2. 

The spectra were made in the same way as in the former cases. The background regions 
were selected from the downstream side of the filaments without other structures. The 
emission is free from thermal photons, according to previous studies, especially in such a 
hard band (Rho et al. 2002). The lower- right panel of Figure 3 shows the background- 
subtracted spectrum, which was accumulated from both filaments. The spectral fittings 
were made with the three models (NEI, power-law, and SRCUT) with an absorption. We 
adopted the spectral index at 1 GHz with the value 0.6 derived by Caswell et al. (1975). 
All of the models well reproduced the data with the best-fit profiles and values given in the 
panel and Table 3. 



4. Discussion 

4.1. Emission Origin of Filaments 

Table 3 summarizes the spectral fittings for the filaments in the four SNRs to thermal 
(NEI) and nonthermal (power-law and SRCUT) models. Because all models give a similar re- 
duced x^, we cannot distinguish which model is the best statistically to reproduce the spectra 
of the filaments. The thermal model fittings for all cases, however, require an unusually high 
temperature and a low abundance for young (ejecta dominated) SNRs. On the other hand, 
the best- fit photon indices in the power-law model are ~2-3, which are similar to hard X- 
rays in other SNRs (e.g., Koyama et al. 1995). It reminds us of the synchrotron emission 
from electrons with a power-law distribution of index of F = 1.5 or softer, suggesting that 
these filaments are acceleration sites of electrons via the DSA mechanism. The fitting of the 
SRCUT model requires a roll-off frequency of about 1 keV, as shown in Table 3. In order to 
check the infiuence of uncertainty of the radio index s, we fitted the spectra with the frozen 
s of 0.5 and 0.9. For the s — 0.5 cases, which is the theoretical minimum indicated by the 
DSA, Vroiioff remains unchanged within a factor of two. On the other hand, the fittings with 
s — 0.9, which implies that the de-acceleration is very efficient, require the change of froiioff 
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in the order of 1. However, the best- fit normahzation in the radio band becomes larger 
than 10% of the total flux from the whole SNR, which is unrealistic for such small regions. 
Moreover, in the Kepler case, the radio map shows that the value of s of the rim is smaller 
than that of the inner region (DeLaney et al. 2002), although it changes from 0.85 to 0.6. 
This result suggests that the outer rims of SNRs are acceleration sites and dc-acceleration 
may not occur yet, in other words, s is around 0.6. For the Tycho case, Katz-Stone et al. 
(2000) indicated that the spectral index of filaments in outer rim is around 0.5, indicating 
that lyroiioff derived in our analysis for the SNR is roughly correct. These results indicates 
that the X-ray photons shown in Figure 1 are emitted by electrons with the highest energy 
in each SNR. 

It can be argued that the thermal spectra in source and background regions may be 
the same. The fine emission from the NEI plasma changes as the ionization time scale (ngt) 
varies. We take the background regions just behind the source regions, and their widths are 
much small (1-5% of SNR radii). The electron density Ue may not change so drastically in 
such thin downstream regions. Since we take the background regions just behind the source 
regions, and their widths are much small (1-5% of SNR radii), the electron density rie docs 
not so drastically change in the downstream regions. Thus, the difference of Uet between the 
background and source regions may not be large enough to clarify the difference between 
their thermal spectra. 

All of the filaments have much smaller scale widths than that derived from the Sedov 
self-similar solution (= R/12), on both the upstream and downstream sides, which is similar 
to the previous results in the SN 1006 case (Bamba et al. 2003b). We should consider the 
infiuence of PSF, especially for the filaments in Cas A and Tycho, where Wu are as small 
as the PSF (see Table 2). The PSF enlarges the apparent Wy, and w^,- Thus, our result 
may be the upper limit of the physical scale width. Together with the fact that our sample 
covers almost historical SNRs with an X-ray shell, it may be suggested that all relative young 
(<2000 years) SNRs accelerate electrons in very thin regions around shock fronts. 

4.2. Comments on Individual SNRs 
4.2.1. Cas A 

We analyzed filaments thinner than that by Vink & Laming (2003), and the former 
has harder spectra than the latter. This fact may imply that the magnetic field around 
our filaments is stronger than that around the filament by Vink & Laming (2003). Based 
on analysis of the filament-like structure in Cas A by Vink & Laming (2003), Berezhko & 
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Volk (2004) estimated the magnetic field around the structure to be ~ 500 fj,G, generated 
by the cosmic rays themselves (Lucek & Bell 2000). Our thinner filaments may have a larger 
magnetic field according to their discussion. 

4-2.2. Kepler's Remnant 

As for Kepler, since there is no report about nonthermal filaments, this is the first 
analysis of them. Cassam-Chenai et al. (2003) reported that south-eastern edge of the SNR 
emits nonthermal X-rays with XMM-Newton. Their spectral parameters are consistent with 
ours. Therefore, the region by Cassam-Chenai et al. (2003) may be a part of our filaments. 

4.2.3. Tycho's Remnant 

Our spectral results are roughly consistent with the results by Hwang et al. (2002), 
although the regions for spectral analysis are slightly different, except for the S abundance 
in the NEI model. The difference may come mainly from selection of the different background 
regions; we chose them not from out of the SNR, but from just at the downstream side of 
the filaments. We could thus extract pure nonthermal emission. 

4.2.4. RCW 86 

We re-confirmed the result of spectral analysis by Rho et al. (2002). A similar dis- 
cussion about the magnetic field in the SN 1006 case by Yamazaki et al. (2004b) may be 
possible, although the filaments are very clumpy and wc may have to consider the curvature 
effect suggested by Berezhko et al. (2003). Together with the shock velocity of 562 km s~^ 
(Ghavamian ct al. 2001) and results of spatial and spectral analysis (sec previous sections), 
we estimated the magnetic field around the filaments is ~ 4-12 /xG (see also Bamba 2004). 
Watanabe et al. (2003) reported on TeV gamma-ray observations of the SW region of this 
remnant using the CANGAROO telescope. According to our results, inverse Compton TeV 
gamma-rays may be detected with the current TeV 7-ray telescopes. The estimated fiux of 
inverse Compton emission is ~ 6 x 10~^^(Sd/10 //G)~^ ergs s~^cm~^. 
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4.3. Correlations between the age and spatial and spectral pEtrameters 

We found thin and nonthermal filaments in four historical SNRs, which may be acceler- 
ation sites of high-energy electrons. Their scale widths and spectral parameters include a lot 
of information about physical parameters inevitable for understanding acceleration on the 
shock, such as the gyro radius of electrons, the magnetic field configuration, and so on, as 
already suggested by Bamba et al. (2003b) in the SN 1006 case. Furthermore, a correlation 
between the age of SNRs and the spatial and spectral parameters might represent their time 
evolution in a single SNR, although we should consider the characteristics of the individ- 
ual source. In this section, we consider some correlation between the age and the physical 
parameters as a rough discussion about the time evolution. 

Figure 4(a) represents the average scale width weighted by the flux of each fllament 
{wu and Wd in the unit of pc) as a function of the age (tage in the unit of year) for five 
historical SNRs; we analyzed four and one was from the previous results of SN 1006 (Bamba 
et al. 2003b). The distances to the SNRs are fixed to previously estimated values (see §3). 
The panel shows that in all SNRs, both Wu and Wd are quite smaller than the SNR radii. 
Moreover, there is a positive correlation between Wd and tage] the scale width becomes larger 
as the SNR becomes older. We fitted them with a power-law function as a tentative model, 
and obtain the relations such as: 

Wd = S.Ot'oi X 10-%,^''-°-^ , (2) 

with the reduced oi 9.68/3. Hereinafter, the errors indicate 90% confidence regions. 
The index is larger than that for the shock width of thermal gas derived from the Sedov 
solution (=2/5), although we reached no definite conclusion because the fittings were rejected 
statistically. 

Figures 4(b) and (c) show the plots of photon indices in the power-law fittings vs. the 
age {tage)0^) and Vroiioff (Hz) of the SRCUT model vs. tage (c), respectively. The spectrum 
seems to grow softer in the very young phase (less than 500 years), whereas the spectral 
parameters remain constant at 1000-2000 years old. Since the error is too large, especially 
in the young phase, we gave up modeling the time evolution of the spectral parameters. 

We searched for another parameter that has a clear correlation with time, and found a 
tentative function, such as 

B = VroUoffWd''^ . (3) 

Figure 4(d) represents the relation between tage and B. It shows monotonous decay, and a 
power-law fitting {B = Ctage") was accepted statistically with best-fit values of 

C = 2.6tiixl0'^ (4) 
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a = -2.96l°:J^ , (5) 

with a reduced of 2.03/3, respectively. This resuh imphes that the function B might 
include some physical quantities that evolve with time. We checked the uncertainty of B 
due to the change of s as already discussed in §4.1, and found that the best- fit value does 
not change although error regions become 10 times larger. Therefore, the uncertainty of s 
cannot change our results. 

Here, we make a scenario in order to explain the time evolution of B. Our assumption 
is that the spatial profile found in our present analysis reflects that of accelerated electrons 
(Bamba et al. 2003b; Yamazaki et al. 2004b). The maximum energy of accelerated elec- 
trons, Emax, is determined by the age of the SNR, or the synchrotron energy loss, such as 
iacc ~ niin{iage, ^loss}, where tacc and tioss are the acceleration and the synchrotron loss time 
scale, respectively (Yamazaki et al. 2004b). Then, we simply obtain the scale width of the 
accelerated electrons in the radial direction, A^^, as A^j ^ v^x minjtago ^loss} ~ ^(i^acci where 
Vd is the downstream fluid velocity. The acceleration time tacc is on the order of K/v^, where 
Vs and K are the shock velocity and the diffusion coefficient, respectively. The quantity K is 
assumed to be proportional to the gyro radius of the accelerated electrons, Vg = Emax/^Bdi 
where is the downstream magnetic field (Drury 1983; Yamazaki et al. 2004b). As a 
result, we find A^; ~ K/vg oc B^r^E^ax^s^^ ■ Because of the projection (curvature) effect, 
observed apparent scale width in the downstream region, wa, differs from A^ (e.g., Berezhko 
& Volk 2004). The ratio Wd/Ad depends on uncertainty of radial profile of accelerated elec- 
tron distribution downstream. It is found that for sufficiently small Ad compared with the 
local curvature radius, the ratio Wd/Ad is about 7, 5, or 1 when wc assume the exponential, 
Gaussian, or the top-hat form of the downstream radial profile, respectively. Hence it may 
be assumed that the ratio w^/A^(~ 3) remains unchanged within a factor of 3, so we find 
oc A^ oc Bd^^EmaxVs^^. The value of Vroiioff is proportional to BdEmax^ (Reynolds & 
Keohane 1999). Therefore, the quantity B is proportional to Bd^Vg^. The time evolution of 
B may represent that of Bd^Vg^. 

Let us consider the time evolution of the magnetic field and the shock velocity as Bd oc 
tage ^ and Vs OC tage The index, m, changes from to 3/5 as SNRs evolve from the 
free-expansion phase to the Sedov phase. Together with the observational fact of eq.(5), the 
index i takes the value of 0.57-1.02. Our result suggests that the magnetic field decreases 
as the SNR ages. This is the first observational implication of the time evolution of the 
magnetic field in the SNR. 

The most interesting case is £ — 0.6, which is allowed when SNRs are in the Sedov phase 
(m = 0.6). In this case, since the density of the thermal plasma is almost constant, the energy 
density of the magnetic field {ub oc Bd'), and the kinetic and thermal energy densities of 
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the shock {uth oc shock temperature and Ukin oc Vs'^), evolves according to the same time 
dependency (index = —1.2). The time evolution of the energy density of accelerated protons 
(up) may also be discussed from the relations by Lucek & Bell (2000) as 

V 

Up = —Ub oc BdVs oc tage^^'"^ , (6) 

Va 

where va = B/^/47Tp and p are the Alfven velocity and the fluid density, suggesting that 
accelerated protons in an SNR evolve with the same energy density evolutions as that of 
the other energy carriers. Since Bamba et al. (2003b) indicated that these energy densities 
are roughly in equipartition with each other, our result implies that they evolve while main- 
taining equipartition. Bell & Lucek (2001) suggested theoretically that the magnetic field 
on the shock evolves according to the same time dependency as the shock velocity, which is 
consistent with our result. 

There remain many unsolved problems that we should consider. The value of Wd depends 
on many parameters that we have ignored, such as the curvature effect, the angle between 
the magnetic field and the shock normal, the degree of turbulence, and so on. Solving these 
problems is beyond our work, which is going to be considered in Yamazaki et al. (2004a). 



5. Summciry 

We have conducted systematic spectral and spatial analysis of filamentary structures in 
historical SNRs for the first time. A summary of our results is as follows: 

1. We discovered very thin filaments on the outer edges of four historical SNRs: Cas A, 
Kepler, Tycho, and RCW 86. The scale widths were measured for the first time; they 
are below or the same as the PSF size of Chandra and 0.01-0.4 pc on the upstream 
and downstream sides of the emission peak, respectively. 

2. Although both thermal and nonthermal spectral models can reproduce the spectra of 
the filaments, the former is unlikely since it requires an unrealistic high temperature 
and low abundances. The photon indices (~2-3) and Vroiioff (~ 10^^-10^*^ Hz) of 
nonthermal models suggest the synchrotron X-ray emission from electrons accelerated 
via DSA, in a similar way to the SN 1006 case. 

3. The value of Wd (see eq.(2)) becomes larger as the SNR ages. 

4. The spectra of filaments might become softer as the SNR ages in the phase of tage ^ 
500 years, whereas it does not change when 1000 < tage ^ 2000 years, although there 
is no definite conclusion because of a lack of statistics. 
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5. The time evolution of a tentative function B — ^roUoff/wd' was analyzed, and it has 
a clear negative correlation with tage with the index of ~ —3. This may imply that 
the downstream magnetic field decreases with the index {€) of 0.6-1.0. When ^ = 0.6, 
energy densities of magnetic field, shock (thermal and kinetic) , and cosmic rays evolve 
while keeping equipartition with each other. 



Our particular thanks are due to M. Hoshino, S. Inutsuka, K. Makishima, and F. Taka- 
hara, for their fruitful discussions and comments. We also thank an anonymous referee for 
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Fig. 1. — Close-up views of the SNRs in the hard X-ray band (written in each image; see 
also text) with Chandra. Grayscales are in logarithmic and coordinates are in J2000. The 
filament regions for spatial and spectral analysis are also shown with solid rectangles. 
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Fig. 2. — Profiles of the filaments of SNRs in the hard X-ray band (Cas A: 5.0-10.0 keV, 
Kepler: 4.0-10.0 keV, Tycho: 2.0-10.0 keV, RCW 86: 2.0-10.0 keV). The best-fit models are 
shown with solid lines. The shock runs from right to left in the all panels. Note that fittings 
were carried out in a part around the shock front in order to avoid any contamination from 
bright thermal emission. 
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Fig. 3. — Spectra of combined-filaments for each SNR (crosses). The best-fit power-law 
models are shown in the solid hnes. The black and gray crosses in the panel for Tycho are 
the spectra accumulated from filaments on the different CCD chips. The lower panel for 
each figure represents the residuals from the best-fit model. 



Table 1. Observation Log 



Target Type ObsID R.A. Dec. Date Exposure (ks) 



Cas A lb 00114 23i^23™40!2 58?47™34!1 2000 Jan. 30-31 50 

Kcplcr ? 00116 17l>30™39!6 -21?30™32!5 2000 .Juri. 30 - ,Iul. 1 49 

Tycho la 00115 10^25™07!0 64'?09™44!7 2000 Sep. 20-21 49 

SN 1006 la 00732 15i'03™51?6 -4ld51?'18?8 2000 Jul. 10-11 68 

SW shell of ROW 86 . II? 01993 UHcrms -62^39^°43?9 2001 Feb. 1-2 92 




Fig. 4. — Time evolutions of parameters, (a) time vs. scale width, (b) photon index of 
power-law fitting, (c) t'roUofr of SRCUT model fitting, and (d) B = I'ronos/wd^. Solid lines 
represent best-fit models for each plot (see text). 



-22- 



Tablc 2. Best-fit scale width of filaments' 



No. A Wu reduced 

[counts arcsec"^] [arcsec] [arcsec] [x^/d.o.f.] 

Cas A 

1 59 (46-72) (<0.93) 1.27 (0.96-1.81) 99.7/61 

2 34 (24-49) (<0.80) (<1.59) 27.3/26 

Kepler 

1 33 (26-41) 1.17 (0.87-1.59) 0.93 (<1.41) 19.3/21 

2 29 (25-34) 1.59 (1.17-2.19) 3.09 (3.46-3.87) 70.1/47 

Tycho 

1 129 (119-139) 1.18 (1.01-1.32) 5.36 (4.47-6.12) 82.7/55 

2 62 (48-70) (<0.80) 1.70 (1.32-3.15) 42.3/32 

3 80 (74-83) (<0.80) 2.38 (2.20-2.54) 23.4/37 

4 150 (146-152) 0.86 (0.80-0.93) 5.53 (5.00-6.14) 48.9/46 

5 63 (57-71) 1.03 (0.90-1.35) 2.47 (1.93-3.15) 17.9/21 

RCW 86 

1 28 (25-30) 2.39 (1.48-3.34) 20.1 (17.3-23.8) 113.4/95 

2 18 (14-23) 1.56 (0.49-4.79) 18.2 (11.8-35.6) 46.8/31 



^Parentheses indicate single parameter 90% confidence regions. 



Table 3. Best-fit parameters of spectral fittings of filaments^ 



NEI power-law SRCUT*' 

No. kT [keV] ^ EM? JVj? xVd-o-f- T JVi? FluxS x^d.o.f. i^.oizo//'^ S^H^^I^ x'/d.o.f. 



Cas A 

1 028 ^^"^ X& 27 28.0/31 2I 27 To 29.2/33 9^ 26 O 29.2/33 

(>2.1) (<l-4) ••• j (2.0-8.1) (15-41) (2.0-3.5) (20-47) (> 4.5) (20-33) (5.1-7.6) q 

2 5.7 ••• 0.030 0.15 17 19.1/27 2.2 20 1.1 22.1/29 28 20 2.1 22.1/29 

(2.7-34) (>0.19) (0.011-0.096) (0.073-3.5) (11-25) (1.9-2.8) (13-31) (> 1.3) (16-24) (1.7-2.4) 

total . 5.4 0.26 3.8 5.6 25 56.4/42 2.1 23 3.1 57.6/44 24 22 6.3 57.6/44 

(3.0-13) (0.019-0.73) (1.6-7.7) (4.1-9.3) (17-35) (1.6-2.7) (15-34) (>1.8) (19-27) (5.4-7.2) 

Kepler 

1 3^5 ~ 8^2 5^5 4^3 43.2/52 Z3 5^8 Tt 42.2/54 5^5 5^ 090 42.8/54 

(2.7-4.9) (< 0.067) (3.4-28) (4.5-6.8) (3.4-5.6) (2.1-2.5) (4.8-7.0) (2.1-17) (4.7-6.0) (0.83-0.98) 

2 3.0 0.13 ••• 7.9 3.7 25.7/30 2.4 5.5 2.2 24.9/32 4.0 4.9 1.4 25.0/32 

(2.5-3.8) (< 0.40) ■ ■ ■ j (6.5-9.7) (3.0-4.5) (2.2-2.6) (4.5-6.7) (1.8-9.8) (4.3-5.5) (1.3-1.5) 

total. 3.2 ••• 11 14 4.0 47.6/53 2.4 5.7 4.1 44.9/55 3.6 5.2 2.8 45.5/55 

(2.6-3.9) (<0.34) (>6.9) (12-17) (3.4-4.7) (2.2-2.6) (4.9-6.6) (2.0-7.9) (4.8-5.7) (2!.6-3.0) 

Tycho CO 

1 Tg 022 3.4 X 10^ 4?7 6^5 102.3/68 3^0 9I 3I 104.8/70 O40 sl 10.70 106.0/70 

(1.6-2.4) (0.062-0.50) (>69) (3.5-6.0) (5.4-7.5) (2.8-3.3) (7.9-11) (0.25-0.70) (7.4-8.8) (0.64-0.77) 

2 1.0 ••• ■•• 2.5 12 9.40/11 3.8 16 0.49 9.05/13 0.11 14 1.4 9.15/13 

(058-6.9) (<1.6) ■■• j (0.27-9.2) (6.6-20) (2.5-5.6) (8.6-27) (0.020-1.3) (11-19) (1.0-2.0) 

3 2.1 0.12 ■•• 1.7 11 40.2/41 2.5 7.6 1.1 43.5/43 1.3 6.9 0.051 43.2/43 

(1.6-2.9) (<0.37) ••• j (0.85-2.4) (5.7-15) (2.2-2.9) (6.0-9.8) (0.52-4.5) (5.9-8.2) (0.045-0.059) 

4 4.3 ■■■ ••• 2.9 5.1 39.5/32 2.1 6.8 4.4 40.9/34 6.5 6.2 0.052 40.6/34 

(3.0-7.0) (<0.38) ••• j (2.3-3.9) (3.4-10) (2.0-2.4) (4.6-9.6) (1.8-46) (5.0-7.7) (0.046-0.059) 

5 1.8 0.05 ••• 2.0 9.8 51.9/47 2.9 9.4 1.1 55.5/49 0.64 8.2 0.12 54.4/49 

(1.4-2.6) (<0.20) ••• j (1.1-3.2) (5.1-14) (2.5-3.3) (6.9-13) (0.26-1.7) (6.7-10) (0.10-0.14) 

total . 2.3 0.07 • • • 12 5.8 182.0/134 2.7 8.0 11 187.7/136 0.83 7.2 0.78 185.2/13f 

(2.0 2.6) (<0.17) ■■• j (11 11) (5.2 6.1) (2.6 2.9) (7.3 9.0) (0..59 1.2) (6.7-7.7) (0.7.'^ 0.8.S) 

i\C\\ «(i 

1 2^5 0.032 ~ 2I 5^7 67.2/58 2I bA lA 68.9/60 2^0 019 68.8/60 

(1.8-3.4) (<0.12) ••• j (1.6-2.8) (3.3-8.3) (2.2-2.8) (4.3-6.9) (0.77-5.9) (4.2-5.7) (0.17-0.22) 

2 1.4 0.27 0.35 0.82 7.0 18.7/25 2.5 2.7 0.32 25.7/27 1.3 2.3 0.051 25.0/27 

(1.0-2.2) (0.019-1.0) (0.20-0.52) (0.49-1.4) (3.2-9.4) (2.1-3.2) (1.1-4.1) (0.29-12) (1.5-3.5) (0.039-0.064) 

total. 2.1 0.11 0.43 3.3 7.0 75.3/85 2.4 4.7 1.8 82.2/87 2.1 4.3 0.22 80.8/87 

(1.6-2.6) (<0.25) (<0.80) (2.1-4.3) (3.1-9.1) (2.2-2.7) (4.2-6.0) (1.0-5.9) (3.7-5.0) (0.20-0.25) 



heses indicate single parameter 90% confidence regions. 

1 index at 1 GHz is assumed according to the previous radio observations (see text), 
ance ratio relative to the solar value (Anders & Grevesse 1989). 

ion time-scale in the unit of 10^" s cm~^, where Ue and t are the electron density and age of the plasma. 

on measure in the unit of lO^^cm"^, with the assumed distance to SNRs according to the previous observations (see text). 

ition column in the unit of lO^^cm"'^, calculated using the cross sections by Morrison & McCammon (1983) with the solar abundances (Anders & Grevesse 1989). 

'ed flux in the 0.5-10.0 kcV band in the unit of 10~^^ ergs cm~^s~^. 

T frequency in the unit of 10^^ Hz. 

jnsity at 1 GHz in the unit of 0. 1 Jy. 

termined. 



